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Abstract 

For the nonlinear wave equation utt—c{u) {c{u)ux) ^ = 0, it is well known that solutions 
can develop singularities in finite time. For an open dense set of initial data, the present 
paper provides a detailed asymptotic description of the solution in a neighborhood of each 
singular point, where \ux\ —t oo. The different structure of conservative and dissipative 
solutions is analyzed. 


1 Introduction 


The nonlinear wave equation 

utt - c{u){c{u)ux)^ = 0 , ( 1 . 1 ) 

provides a mathematical model for the behavior of nematic liquid crystals. Solutions have 
been studied by several authors [1, 2, 3, 4, 7, 10, 12, 13]. We recall that, even for smooth 
initial data 

u{x,0) = uo{x), ut{x,0) = ui{x), (1.2) 

regularity can be lost in finite time. More precisely, the norm of the solution u(-, t) remains 
bounded, hence u is always Holder continuous, but the norm of the gradient ||rta;(-,t)||L°° can 
blow up in finite time. 

The paper [4] introduced a nonlinear transformation of variables that reduces (1.1) to a semi- 
linear system. In essence, it was shown that the quantities 

w = 2 avctan{ut + c{u)ux), z = 2 avctan[ut — c{u)ux), 

satisfy a first order semilinear system of equations, w.r.t. new independent variables X, Y 
constant along characteristics. Going back to the original variables x, t, u, one obtains a 
global solution of the wave equation (1.1). 
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Based on this representation and using ideas from [5, 6, 8, 9], in [1] it was recently proved 
that, for generic initial data, the conservative solution is smooth outside a hnite number of 
points and curves in the t-x plane. Moreover, conditions were identified which guarantee the 
structural stability of the set of singularities. Namely, when these generic conditions hold, the 
topological structure of the singular set is not affected by a small perturbation of the initial 
data. 

Aim of the present paper is to derive a detailed asymptotic description of these structurally 
stable solutions, in a neighborhood of each singular point. This is achieved both for conser¬ 
vative and for dissipative solutions of (1.1). We recall that conservative solutions satisfy an 
additional conservation law for the energy, so that the total energy 

£{t) = ^ j[uf + c^{u)ul]dx 

coincides with a constant for a.e. time t. On the other hand, for dissipative solutions the total 
energy is a monotone decreasing function of time. A representation of dissipative solutions in 
terms of a suitable semilinear system in characteristic coordinates can be found in [3]. 

The remainder of this paper is organized as follows. In Section 2 we review the variable 
transformations introduced in [4] and the conditions for structural stability derived in [1]. 
Section 3 is concerned with conservative solutions. In this case, for smooth initial data the 
map 

{X,Y) ^ (x,t,u,w,z)(X,Y) (1.3) 

remains globally smooth, on the entire X-Y plane. To recover the singularities of the solution 
u{x, t) of (1.1), it suffices to study the Taylor approximation of (1.3) at points where rc = tt or 
z = Tx. In Section 4 we perform a similar analysis in the case of dissipative solutions. This case 
is technically more difficult, because the corresponding semilinear system has discontinuous 
source terms. 

We remark that, for conservative solutions, a general uniqueness theorem has been recently 
established in [2]. On the other hand, for dissipative solutions no general result on uniqueness 
or continuous dependence is yet known. Whether structurally stable dissipative solutions are 
generic, arising from an open dense set of initial data, is also an open problem. 


2 Review of the equations 


Throughout the following, on the wave speed c we assume 


(A) The map c : M i—> is smooth and uniformly positive. The quotient c'{u)/c{u) is 

uniformly bounded. Moreover, the following generic condition is satisfied: 


c'{u) = 0 


c"{u) / 0. 


( 2 . 1 ) 


Because of (2.1), the derivative c'(u) can vanish only at isolated points. 

In (1.2) we consider initial data {uo,ui) in the product space ff^(M) x L^(R). It is convenient 
to introduce the variables 

( R = ut + c{u)ux , 2 ) 

\ S = ut- c{u)ux , 
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In a smooth solution, B? and satisfy the balance laws 

' {R\-{cR% = f^{R^S-RS^), 

(2.3) 

(52)i + (c52)^ = ^^{S^R-SR^). 

As a consequence, the energy is conserved: 

1 , 2 2 2\ + 5*^ 

E + = - - -, (2.4) 

One can think of R? and 5^ as the energy of backward and forward moving waves, respectively. 
Notice that these are not separately conserved. Indeed, by (2.3) energy can be exchanged 
between forward and backward waves. 

A major difficulty in the analysis of (1.1) is the possible breakdown of regularity of solutions. 
Indeed, even for smooth initial data, the quantities Ux,ut can blow up in finite time. To deal 
with possibly unbounded values of R, S, following [4] we introduce a new set of dependent 
variables: 

w = 2arctani?, 2 : = 2arctan5. (2-5) 



Figure 1: The backward and forward characteristic through the point {x,t). 


To reduce the equation (1.1) to a semilinear one, it is convenient to perform a further change 
of independent variables (Fig. 1). Consider the equations for the forward and backward 
characteristics: 

x'^ = c{u ), x~ = —c{u). (2.6) 

The characteristics passing through the point (x, t) will be denoted by 

s x'^{s,x,t), s^x~{s,x,t), 


respectively. As coordinates {X,Y) of a point {x,t) we shall use the quantities 


X = X {0,x,t), Y = —x~^{0,x,t). 


For future use, we now introduce the further variables 


P = 


l + i?2 


q = 


1 + 52 


(2.7) 


( 2 . 8 ) 


3 




Starting with the nonlinear equation (1.1), using X,Y as independent variables one obtains a 
semilinear hyperbolic system with smooth coefficients for the variables u, w, z, p, q, x, t, namely 


ux = 

SlllR? ^ 

4c P 


sin 2 

Uy = 

4c y 


(2.9) 


wy = ^ (cos 2 : — cos w) q , 


= 


C 

Sc‘ 


- (cos w — cos z) p , 
(sin z — sin w) pq , 


PY = ^ 

Qx = ^ (sin w — sin z) pq , 


_ (l+cos«))p 

4 

. (l-\-cos z)q 

k XY = 


_ (H-COS117)p 

~ 4c 

_ (H-cos^)^ 

~ 4^ 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


See [4] for detailed computations. Boundary data can be assigned on the line 70 = {{X,Y ); X+ 
y = 0}, by setting 


u[s, —s) = u[s), 

x{s, —s) = x(s), 

t{s, -s) = t{s ), 


w{s, —s) = w{s ), 

z{s, -s) = z{s ), 


p(s, -s) = p{s), 

q{s, -s) =q{s), 


(2.14) 


for suitable smooth functions u, ,x,t,w, z, p, q. 


Remark 1. The above system is clearly invariant w.r.t. the addition of an integer multiple of 
27r to the variables w,z. Taking advantage of this property, in the following we shall regard w, z 
as points in the quotient manifold T = M/27rZ. As a consequence, we have the implications 


W ^ TT 
Z 7^ TT 


COS w > — 1 , 

cos z > — 1 . 


(2.15) 


Remark 2. The system (2.9)-(2.13) is over deter mined. Indeed, the functions u,x,t can be 
computed by using either one of the equations in (2.9), (2.13), ( 2 . 12 ), respectively. As shown 
in [ 1 ], in order that all the above equations be simultaneously satisfied along the line 70 one 
needs the additional compatibility conditions 



sinr(;(s) 

4c(ll(s)) 


p{s) 


sin z(s) 
4c(ll(s)) 


q{s), 


(2.16) 
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(1 + cos q{s))p{.s) + (1 + cos z{s))q{s) 

(2.17) 

4 

d , 

(1 + cosw{s))p{s) — (1 + cos z{s))q{s) 

(2.18) 

4c(ll(s)) 


In turn, if (2.16)-(2.18) hold along 70 , then a unique solution to the system (2.9)~(2.13) can 
be constructed, on the entire X-Y plane. 


Given initial data (t(o,rti) in (1.2), we assign boundary data (2.14) on the line 70 , by setting 


f p(x) = 1 + R^(x,0), 

\ q(x) = 1 + S^(x,0). 

(2.19) 


U[X) = Uo{X) , 

t(x) = 0 , 
x(x) = X , 


w(x) = 2 arctan i?(x, 0 ), 
z(x) = 2 arctan 5(x, 0 ), 


We recall that, at time t = 0, by (1.2) one has 

R(x,0) = (ut + c{u)ux)ix,0) = ui{x) + c{uoix))uo,xix), 

S{x,0) = {ut - c{u)ux)ix,0) = Ui{x) - c{uo{x))uo,xix). 

As proved in [ 1 ], for any choice of uo,ui in (2.19) the compatibility conditions (2.16)-(2.18) 
are automatically satisfied. 

The following theorems summarize the main results on conservative solutions, proved in [4, 1, 
2 ]. As before, U denotes the product space in ( 2 . 21 ). 


Theorem 1. Let the wave speed c(-) satisfy the assumptions (A). 

Given initial data (tto,rti) G H^{W) x L^(R), there exists a unique solution 

{X, Y) {u, w, z,p, q, X, t){X,Y) to the system (2.9)-(2.13) with boundary data (2.14), (2.19) 

assigned along the line 70 . Moreover, the set 

Graph(u) = | (x(A, T), t(A, T), u(A, T)) ; (A,y)GR2| (2.20) 

is the graph of the unique conservative solution u = u{x,t) of the Cauchy problem (1.1)-(1.2). 


Theorem 2. Let the assumptions (A) be satisfied and let T > 0 be given. Then there exists 
an open dense set 


V C U = (c^(R) n X (^C2(R) nL2(R)) 


( 2 . 21 ) 


such that the following holds. 

For every initial data {uo,ui) G V, the corresponding solution {u,w, z,p,q, z,t) of (2.9)-(2.13) 
with boundary data (2.14), (2.19) has level sets {w = tt} and {z = tt} in generic position. 
More precisely, none of the values 


{w,wx,wxx) = (vr, 0 , 0 ), 
{z,zy,zyy) = (vr, 0 , 0 ). 


( 2 . 22 ) 
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(w,z,wx) = 
(w,z,zy) = 
(w,wx,c'(u)) 
(z,ZY,c'(u)) 
is ever attained, at any point {X,Y) for which 

{x{XX), t{X,Y)) 


(7r,7r,0), 


(2.23) 

(7r,7r,0), 


= (7r,0,0), 

(2.24) 

N 

0 


G R X [0, T] . 

(2.25) 


The singularities of the solution u in the x-t plane correspond to the image of the level sets 
{w = tt} and {z = n} w.r.t. the map 

A: (X,y) ^ (x(x,y), t(x,y)). ( 2 . 26 ) 

If none of the values in (2.22)-(2.24) is ever attained, by the implicit function theorem the 
above level sets are the union of a locally finite family of curves in the X-Y plane. In turn, 
restricted to the domain M x [0,T], the singularities of u are located along finitely many 
curves in the x-t plane. 



Figure 2: Two level sets {w = tt} and {z = tt}, in a generic conservative solntion of (2.9)-(2.11). 
Here P is a singular point of Type 1, while Pi, P2, P3 are points of Type 2, and Qi, Q2 are points 
of Type 3. Notice that at Pi, structural stability requires that the function Y{X) implicitly defined 
by w{X,Y{X)) = TT has strictly positive second derivative. At the points Qi, Q2, by (2.10) one has 
wy = zx = 0- Hence the two curves {ru = tt} and {z = tt} have a perpendicular intersection. 


3 Singularities of conservative solutions 


For smooth data uo,ui G the solution {X,Y) i-)- {x,t,u,w, z,p,q){X,Y) of the semi- 

linear system (2.9)-(2.I3), with initial data as in (2.14), (2.19), remains smooth on the entire 
X-Y plane. Yet, the solution u = u{x,t) of (1.1) can have singularities because the coordinate 
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Figure 3: The images of the level sets {w = tt} and {z = tt} in Fig. 2, under the map A : (X, Y) >->■ 
{x{X, Y), t{X, F)). In the x-t plane, these represents the curves where u = u{x, t) is not differentiable. 
A generic solution of (1.1) with smooth initial data remains smooth outside finitely many singular 
points and finitely many singular curves, where Ux —>■ ±oo. Here pi,P 2 ,P 3 are singular points where 
two new singular curves originate, or two singular curves merge and disappear. At the points qi,q 2 a 
forward and a backward singular curve cross each other. 


change A : {X,Y) i—>■ {x,t) is not smoothly invertible. By (2.13)-(2.12), its Jacobian matrix is 
computed by 


DA 


f XX 

XY \ 


/ (l+costt;)p 

I 4 

(l+cos2:)g \ 

V 

ty J 


[ (H-COStl7)p 

(l+cos 2 :)g / 




\ 4c('u) 

4c('u) / 


(3.1) 


We recall that p, q remain uniformly positive and uniformly bounded on compact subsets of 
the X-Y plane. By Remark 1, at a point (XojYo) where w ^ n and z ^ tt, this matrix 
is invertible, having a strictly positive determinant. The function u = u{x,t) considered at 
(2.20) is thus smooth on a neighborhood of the point 


{xo,to) = (x(Xo,Fo), 

To study the set of points x-t plane where u is singular, we thus need to look at points where 
either w = n or z = n. 


If the generic conditions (2.22)~(2.24) are satisfied, then we have the implications 

-(cos z + l)q ^ 0 , 

-(cos w + l)p yf 0 . 
oc“ yu j 

Therefore, by the implicit function theorem, the level sets 


and wx = 0 = 

c'(u) 

^ = 8cy») 

and zy = 0 = 

c'(u) 

^ Sc^ (u) 


A {(X,y); w{X,Y) =7r}, S^X{{X,Y)-, z{X,Y) = 7r}, (3.2) 


are the union of a locally finite family of smooth curves. The singularities of u in the x-t plane 
are contained in the images of S'^ and under the map (2.26). Relying on Theorem 2, we 
shall distinguish three types of singular points P = (Xo,Fo)- 
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( 1 ) Points where w = tt but wx ^ 0 and z / vr (or else, where z = tx but zy 7 ^ 0 and 
w ^ tt). 

(2) Points where w = tx and wx = 0, but wxx / 0 (or else: z = tx and zy = 0, but 

/ 0 ). 

( 3 ) Points where w = tx and z = tx. 

Points of Type 1 form a locally finite family of curves in the X-Y plane (Fig. 2). Their 
images A(P) yield a family of characteristic curves in the x-t plane where the solution u = 
u{x,t) is singular (i.e., not differentiable). 

Points of Type 2 are isolated. Their images in the x-t plane are points where two singular 
curves initiate or terminate (Fig. 3). 

Points of Type 3 are those where two curves {rc = vr} and {z = tt} intersect. Their image in 
the x-t plane are points where two singular curves cross, with speeds ±c{u). 


Our main result provides a detailed description of the solution u = u{x, t) in a neighborhood of 
each one of these singular points. For simplicity, we shall assume that the initial data (rto,rti) 
in (1.2) are smooth, so we shall not need to count how many derivatives are actually used to 
derive the Taylor approximations. 


Theorem 3. Let the assumptions (A) hold, and eonsider generic initial data {uo,ui) G V 
as in (2.21), with uo,ui G C'°°(M). Call {u,w, z,p,q,x,t) the corresponding solution of the 
semilinear system (2.9)-(2.13) and letu = u{x,t) he the solution to the original equation (1.1). 
Consider a singular point P = (Ao,To) where w = tx, and set (xo,to) = (^^(Ao, Fb); ^(Ao, To)). 


(i) If P is a point of Type 1, along a curve where w = tx, then there exist eonstants a ^ 0 
and bi , 62 such that 

I 2/3 


u{x,t) = u{xo,to)-a- c{uo){t - to) {x - Xo) 


+bi ■ {x - Xo) + 62 • (t - to) + 0 ( 1 ) • (^\t - to| + \x- xol) 


4/3 


(3.3) 


(a) If P is a point of Type 2, where w = tx, wx = 0 , and wxx > 0 , then in the x-t plane 
this eorresponds to a point (xo,to) where two new singular eurves 7~,7^ originate. In 
this ease, there exists a constant a ^ 0 such that 

X 1 \ 

u{x,t) = u{xo,to) -I a ■ c{uo){t-to)-I {x - Xo) + 0 ( 1 ) • ( |t - to| + |a: - xo| j 

(3.4) 

(Hi) If P is a point of Type 3, where w = z = tx, then in the x-t plane this corresponds to a 
point {xo,to) where two singular eurves 7,7 cross each other. In this case, there exist 
constants oi 7^ 0 and 02 7^ 0 such that 

1 1 

u{x,t) = tt(xo,to) + Oi • c(tio)(t - to) + (a^ - a^o) 


+02 


c(r(o)(t - to) - (a: - Xo) ^ + 0(1) • (^|t - to| + |x - xo|^. 


(3.5) 



Throughout the following, given a point P = {Xq,Yq) in the X-Y plane where w = tt, we 
denote by {uo,wo, zo,po, qo, xo,to) the values of {u,w, z,p,q,x,t) at {Xq^Yq). The three parts 
of Theorem 3 will be proved separately. 


3.1 Singular curves. 


Let P = (Xo,To) be a point of Type 1, where 

wq = tt, zo / TT, wx{Xo,Yo) 7^ 0. (3.6) 


By the implicit function theorem, the level set where re = tt is locally the graph of a smooth 
function X = ‘h(T), with <h(To) = ^o- We claim that, in a neighborhood of the point 
(xo,to) = the image A(5“') is a smooth curve in the x-t plane, say 

7 = {(a;T); X = (j){t)}. ( 3 . 7 ) 

Indeed, the curve 7 is the image of the smooth curve {X = ^>(y)} under the smooth, one-to- 
one map 

F ^ (x(<h(F),F), t(<h(F),y)). 

For future record, we compute the first two derivatives of (/> at t = to- Differentiating the 
identity rc(<h(F),y) = tt one obtains 

wx^' + VJy = 0 , 


Wxx ■ + WYY + = 0. 

By (2.12)-(2.13), at the point (XojFo) we have 

(1 + coszo)q'o (l + coszo)go 


^{xmY),Y), mY),Y)) = (- 
Observing that 

</.'(t(d>(y),y)) = 

at t = to we have 


4 ’ 4c(uo) 

xx(^(y), y) • ^'(y) + xy ($(y), y) 

tx(^(Y),Y)-^'(Y) + tY(^(Y),Y) ’ 
</>'(*o) = - c(uo ). 


^ ( 0 , 0 ). 


In a similar way we find 

XYY(Xo,Yo)tY(Xo,Yo) — tYY(Xo,Yo)xY(Xo,Yo) c'(uo)smzo 


rito) = 


tUXo,Yo) 


1 + cos Zo 


Next, by (2.9) one has 

ux{Xo,Yo) = 0, uy{Xo,Yo) = go = ai. (3.8) 

Differentiating the first equation in (2.9) w.r.t. X and using (2.10)-(2.11) we obtain 


uxx 


cos w sin w , sin w n sin w 

TTX^^P ~ A 2! ■ ITTYP + ITTYP^ ’ 

4c(u) 4c"'(ii) 4c(u) 4 c(m) 
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uxx{Xo,yo) = = 0^2 7^ 0, (3.9) 

uxxx{Xo-,yQ) = - -^^^^^{wxx{Xo-,yo) Pq+ ‘^wx{Xo,Yq)px{Xq,Yq)^ = as, (3.10) 

«y(X„.r„) = -^.|(g_(cosi 04 . (3,11) 

This yields the local Taylor approximation 

u{X,Y) = uo + ai{y- To) + y (^ " ^o)' + Y + «4 (^ - X^W - Tq) 


+0(1) • (|X - Xol^ + |y - Top + Y - Top |T - Tol). 

(3.12) 

Using (2.13), we perform an entirely similar computation for the function t in a neighborhood 
of (To, To). 

tx(To,To) = 0, ty(To,To) = = /3i > 0, (3.13) 


sin w 1 + cos w, 1 + cos w 

txx = - —rxVJxp - TYr~r'^ W ^- A ( \ P^ ’ 

4c(u) 4c^(u) 4c(u) 

ixx(To,To) = txy(To,To) = 0, 

o (V w\{Xo,Yo) . o , f. 

txxx(Xo,Yo) = —— P^ = /^s + 0 . 

This yields the Taylor approximation 


(3.14) 

(3.15) 


t{X,Y) = to + (3,{y -yo) + ^{X -Xof 

b 

+0(1) • (|T - Top + |T - Top + |X - Xop |T - To|). 


(3.16) 


Finally, for the function x, using (2.12) we find 

xx(To,To) = 0, xy{Xo,Yo) = -- - qo = - 71 < 0 , (3.17) 

sin w ■ wx 1 + cos w 
XXX = -^-pH-^- Px , 


xxx{Xo,Yq) = 0, xxy{Xo,Yo) = 0, (3.18) 

xxxx(Xo,Yo) = po = 73 > 0. (3.19) 

This yields the Taylor approximation 

x(X,T) = xo - 71 (T - To) + 73 (T - To )3 

(3.20) 

+ 0 ( 1 ) • (|T - Top + |T - Top + |T - Top |T - To|). 
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Observing that the above Taylor coefficients satisfy 

7i = c(iio)A, 73 = c { uo )/ 33 , 

from (3.16) and (3.20) we deduce 
(x - Xo) - c{uo){t - to) 

= - 271 {Y - To) + 0(1) • (l^ - + 1^ - + 1^ - - Yo\), 

(x - Xo) + c{uo){t - to) 

= 273 (X - Xo)3 + 0(1). + |y - Tol' + |y - |y - yol) ■ 


(3.21) 


(3.22) 


Next, using (3.16) and (3.20) we obtain an approximation for X,Y in terms of x,t, namely 
1 + cos Zo 


%{Y-Yo) = c{uo){t-to)-{x-xo)+0{l)-{\X-Xo\^+\Y-Yo\^ + \X-Xo\‘^ l^-^o|). 
■po{X-Xo)^ = c{uo){t-to)+{x-xo)+0{iy(\X-Xo\^+\Y-Yo\^+\X-Xo\^ l^->"o|). 


wUXo,Yo: 

12 

Inserting the two above expressions into (3.12), we finally obtain 

9po 


u{t,x) = u{to,Xo)- 


+ 


32wxiXo,Yo 
sinzo 


r -|2/3 

c{uo)it-to) + (x - Xo) 


2c(uo)(l + coszo) 
wxx{Xo, yo) 


c{uo){t - to) - (x - Xo) 


(3.23) 


2ciuo)wl{Xo,Yo) L 

+ 0(1) • (^\t - to| + \x- xol) 
This proves (3.3), with 


c{uo){t - to) + (x - Xo) 
4/3 


a = 


9po 


1/3 


/ 0 . 


^32it;x(-^o,yo), 

The coefficients 61,62 can also be easily computed from (3.23). 


(3.24) 


Remark 3. By (3.3), the solution u is only Holder continuous of exponent 2/3 near the 
singular curve 7 in (3.7). In particular, the Cauchy problem 

x{t) = - c{u{t,x{t))), x{to) = (pito), 

has a solution t i-7> x{t) which crosses 7 at the point (xo,to)- Calling S{t) = x{t) — 4>{t), to 
leading order one has 

5 = c'{uo) ■ . 

Hence, for t ^ to we have 

T'(uo) a'''^ 


6{t) 


{t - to) • 


(3.25) 


The singular curve 7 is thus an envelope of characteristic curves, which cross it tangentially 
(see Fig. 4). 
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Figure 4: Left: a singular curve where w = tt, in the X-Y plane. Vertical lines where X = constant 
correspond to characteristic curves of the wave equation (1.1), where x = —c{u). Right: the images 
of these curves in the x-t plane, under the map A at (2.26). The singular curve 7 is an envelope of 
characteristic curves, which cross it tangentially. 

3.2 Points where two singular curves originate or terminate. 

Let P = (Xo,yo) be a point of Type 2, where 

wq = tt, zo / vr, wx{Xo,Yo) = 0, wxx{Xo,Yo) / 0. (3.26) 


Recalling (2.15), by (3.26) we have 


wy { Xo , Yo ) = + cos 2:0)90 7^ 0. 


By (2.9), at the point (Xo,To) we have 


ux = uxx = 0 , 


UXY = - zo)poqo , 


sin zq 

wxxiXo, To) 

uxxx = - 77 —T- Po 

4c(uo) 


In this case, the Taylor approximation for u near the point (Xo,Lb) takes the form 

(YV\ , tv V\ Wxx{Xq,Yq) 3 

= ““ + 477) * ■ 24c(u,) 


+0(1) - MX - Xol" + |y - Top + |T - Xol |T - To 


(3.27) 


Computing the partial derivatives of x(X,Y) and t{X,Y) at the point (Xo,To), by (2.10) and 
(3.26) we find 


XX = XXX = xxxx = xxxxx = XXY = XXXY = 0, 


(3.28) 


3rc^^(Xo,To) 1 + coszo , on^ 

xxxxxx = -^- Po / 0, xy = -^- go / 0. (3.29) 


ix = txx = txxx = txxxx = txY = txxY = 0 , 


(3.30) 
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txxxxx = 


^o) 


y\. y\. y\. a / \ X^u 

4c(iio) 

This yields the Taylor approximations 


Po 7^ 0 , 


1 + coszo 


(3.31) 


(v\^\ 1 + coszo . ^ I 3wi:x(^0)^o) ^ ^5 

= xo - - - qo(X - To) H-- Po[X - Xof 


+0(1)- (|X-Xo|6 + |T-To|2). 


+0(1)- (|X-Xo|® + |T-To|2), 

Combining (3.32) with (3.33) we obtain 

(X-Xo)^ = 2 y V i -[c(uo)(t-to) + (a;-a;o)] 

3^xx(^o,>o)Po 


(3.32) 


(3.33) 


+ 0 ( 1 ). (|X-Xo |6 + |T-To| 2 ), 


(3.34) 


T-Tn = 


(1 + coszo) qo 


■ [c{uo)it - to) - (x - xo)] + 0(1) • (|X - Xo|6 + |T - Top). 


Inserting (3.34)-(3.35) into (3.27) we eventually obtain 

/ , , , 1 / 803p2 \V5 

u{t,x) = «(to,Xo)-7T7-7-T' T—-77^-i7T ' 


24c(uo) Vrexx(Xo,To) 


c{uo){t - to) + (x - Xo) 


(3.35) 


(3.36) 


+0(1) • (^\t - tol + |x - xol) 


This proves (3.4). 

It remains to show that two singular curves originate or terminate at the point (xq, to)- To fix 
the ideas, assume that 

. u;xx(Xo,To) 

" = - 2 w(x„.r„) > “■ 

By the implicit function theorem, the curve where w = tt can be approximated as 

T-To = k(X-Xo)2 + 0(1)-|X-XoP. (3.38) 


On the other hand, by (3.33) we have 


T-To = «(t-to)-/3(^-^o)^ + 0(l)-(|X-Xop + |t-toP + |^-^o||t-to|), (3.39) 


Muq) . 

(l + coszo)g'o 

Combining (3.38) with (3.39) we obtain 


5! (1 + coszo)g'o 


k{X-Xo)^ = ait - to) + 0(1) ■ \X - Xof. 

Therefore, as shown in Fig. 5 in a neighborhood of (Xo,To) the following holds: 


(3.40) 
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• The two curves {t{X, Y) = tg} and {w{X, Y) = n} intersect exactly at the point (Xq, Tq)- 

• When r < to, the curves {t{X,Y) = r} and {w{X,Y) = tt} have no intersection. 

• When r > to, the curves {t{X,Y) = r} and {w{X,Y) = vr} have two intersections, at 
points Pi = {Xi,Yi) and P 2 = {X 2 ,Y 2 ) with 


Xi — Xo — — J— {t — to) + 0{1) ■ {r — to). 


X 2 -X 0 = +J-{T-to) + 0{l)-{r-to). 


(3.41) 


( Yi-Yo = a{T-to) + 0{l)-{T-to)^/^ 
\ y2-Xo = a{T-to)+ 0 { 1 ) ■ {t- to)^^^. 


w= K 



Figure 5; Left: the equation w(X,Y) = tt implicitly defines a function Y{X) with a strict local 
minimum at Xq. Under generic conditions, y"(Xo) > 0. The dotted curves where t{X,Y) = r have 0, 
1, or 2 intersections respectively, if r < to, r = to, or r > to- Right: the image of the curve {re = tt} 
under the map A in (2.26) consists of two singular curves 7“, 7^ starting at the point po = (xo,to)- 
For T > to, the distance between these two curves is 7 '''(t) — 7 “(t) = 0(1) • (r — to)^'^^. 

For t > to, the solution u = u{t,x) is thus singular along two curves 7~,7''' in the x-t plane 
(see Fig. 5, right). Our next goal is to derive an asymptotic description of these curves in a 
neighborhood of the point (xo,to), namely 

f l~{t) = Xo - c{uo){t-to) + a{t-to)^ - P{t-to)^/"^ + 0{1) ■ {t-to)^, 

{ . (3.43) 

[ 7 +(t) = Xo - c{uo){t-to) + a{t-to)^ + P{t-to)^^"^ + 0 { 1 ) ■ {t-to)^, 

for suitable constants a, /S. 

To prove (3.43), we need to compute more accurate Taylor approximations for t and x near 
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the point (Xojlo)- 


+0(1) • (|X - Xo|6 + |y - FoP +1^ - ^op - ^op), 


(3.44) 


x{X,Y) = xo- 


= - F„) + ky - y,>f + - X„)^ 




■po(^-^o)(^-ho)' 


(3.45) 


+0(1) • {\X - Xol® + |y - Fol" + 1^ - 1^ - ^op). 


The constants a and b are here given by 

zysinz 1 + cosz / c'(«)sinz 

^ = ~ 8ciu) ^ 8ciu) ” 4c2(n) ) ’ 

, 1 1 / 
b = -zyqsmz —-{1 + cos zjqY ■ 

8 8 

where the right hand sides are evaluated at the point {Xq, Iq)- 

For a fixed r > to, let Pi = (Xi,yi) and P 2 = (^ 2 )^ 2 ) be the two points where the curves 
{t{X,Y) = r} and {w{X,Y) = tt} intersect. Let x = 7 “(t) and x = 7 '''(t) describe the 
corresponding points in the x-t plane (see Fig. 5). 

At the intersection point Pi = {Xi,Yi), using (3.41) and (3.42) we obtain 
x(Ai,yi) - Xo + c(uo)(r - to) 

= (a + b)(Y, - y„)= + - Jf„)= + !!iA+»)p„(x. - x„)(n - y„)^ 


+0(1) • (|Ai - Aol^' + - hbl" + |Ai - Aol" - hoP) 


2 f3wj,^{Xo,Yo) , wUXo,Yoy 


= «^(<. + (,)(r - - ‘bY + 0 ( 1 ) ■ (r - <o)’. 

(3.46) 

This yields the equation for 7 “ in (3.43), with suitable coefficients d,l3. An entirely similar 
argument yields the equation for 7 +. In particular, the distance between these two singular 


curves is 


7 '^(i )-7 {t) = 2y{t-tof/^ + 0{l)-\t-tof. 


(3.47) 


15 



3.3 Points where two singular curves cross. 


We now consider a point P = (Xq, lo) where w = z = tt. 

For a generic solution, satisfying the conclusion of Theorem 2, this implies 

wx{Xo,Yo) / 0, zy{Xo,Yo) / 0. (3.48) 

On the other hand, (2.10) yields 

w(Xo,yb) = zx{Xq,Yq) = 0. 

By (2.9) and (2.16), we know that 


ux{Xq,Yq) = Uy{Xq,Yq) = uxy{Xq,Yq) = 0. 
Hence, in a neighborhood of (Xo,To) the function u can be approximated by 


u(X,Y) = uo- 


Wx{Xq,Yo) 2 Zy(Xq,Yo) 


8c{uo) 


'-poiX-Xoy-- 


8c{uo) 


qo{Y-Yoy + 0{l)-{\X-Xo\ + \Y-Yo\y. 

(3.49) 


In addition, by (2.12)-(2.13) we have 

wl{Xo,Yo: 


t{X,Y) = to + 


24c(uo) 


■m(x - x„f + - Y„f 


24c(iio) 


(3.50) 


+0(1)- |x-Xo| + |y-yo 


I(.Y, Y) = J.„(X -x„f- - Y,f 


24 


24 


+0(1)- |x-Xo| + |y-yol . 


(3.51) 


Using (3.50)-(3.51) in (3.49) we eventually obtain 


u{t,x) =u{to,xo) 


1 

8c(uo) 


( 144 Po 

Vu;x(Xo,yo)y 


c{uo){t - to) + {x- Xo)]^^^ 


^ (|t - tol + |x - a:„|). 

(3.52) 

This proves (3.5). □ 


4 Dissipative solutions 

In this last section we assume c'{u) > 0 and study the structure of a dissipative solution in a 
neighborhood of a point where a new singularity appears. We recall that dissipative solutions 
can be characterized by the property that R, S in (2.2) are bounded below, on any compact 
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subset of the domain {(t, x); t > 0, x £ M}. As proved in [3], dissipative solutions can be 
constructed by the same transformation of variables as in (2.5), (2.7), and (2.8). However, the 
equations (2.10)-(2.11) should now be replaced by 


Wy 

= 9 

c'{u) 

8c-‘{u) 

■ (cos Z 

— cos w) q 

< 







= 9- 

c'(u) 

8c^{u) 

(cos w 

— cos z' 

)P, 

PY 

< 

= 9- 

c'{u) 
8c‘‘ (u) 

[ sin z - 

- sin re] 

pq 


= 9- 

c'(u) 

8c^(u) 

[ sin w 

— sin z] 

pq 


where 

^ _ f 1 if max{rt;, z} < tt , 

( 0 if max{rt;, z} > vr . 

Notice that, by setting 9 = 1, one would again recover the conservative solutions. 


(4.1) 


(4.2) 


(4.3) 


It is interesting to compare a conservative and a dissipative solution, with the same initial data. 
Consider a point P = {Xq, Yq) of Type 2, where two new singular curves 'y~ originate, in 
the conservative solution. To fix the ideas, assume that the singularity occurs in backward 
moving waves, so that R —)■ +oo but S remains bounded. Moreover, let the conditions (3.26) 
and (3.37) hold. 

Up to the time to = t{Xo,Yo) where the singularity appears, the conservative and the dissi¬ 
pative solution coincide. For t > to, they still coincide outside the domain 

n = {{x,t); t>to, 'y~{t) < X <j{t)}, (4.4) 

where 7 is the forward characteristic through the point (xo,to)- Figure 6 shows the positions 
of these singularities in the X-Y plane and in the x-t plane. Figure 7 illustrates the difference 
in the profiles of the two solutions for t > to- Our results can be summarized as follows. 


Theorem 4. In the above setting, the conservative solution •) has two strong singu¬ 
larities at X = 7“(t) and x = ^~^{t), where —>■ 00, and is smooth at all other points. 

On the other hand, the dissipative solution •) has a strong singularity at x = 7“(t), 

where 00 , and a weak singularity along the forward characteristic x = 7(t), where 

y(hss continuous but the second derivative does not exist. 

The difference between these two solutions can be estimated as 

|h“"^(U-)-«"*^*(C-)llco(R) = O{l)-{t-to). (4.5) 


Proof. 1. To fix the ideas, assume that at the point P = {Xo,Yo) where the singularity is 
formed one has 

wxx < 0, wy > 0, c\u) > 0. 

In the X-Y coordinates, for smooth initial data the components {x,t,u,w, z,p,q) of the con¬ 
servative solution remain globally smooth. On the other hand, for a dissipative solution by 
(4.I)-(4.2) we only know that these components are Lipschitz continuous. 
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Figure 6: The positions of the singularities in the X-Y plane and in the x-t plane. This refers to a 
point where a new singularity is formed, in the first family (i.e., for backward moving waves). Above: 
a conservative solution. Below: a dissipative solution. Notice that the entire region between the curves 
a~ and a'^ is mapped onto the single curve 7“. Indeed, horizontal segments in the X-V plane are 
mapped into a single point. In the x-t plane, the two solutions differ only on the set fi, bounded by 
the characteristic curves 7“ (the image of both cr“ and cr**) and 7 (the image of the line a). 



T(t) 


Figure 7 : Comparing a conservative and a dissipative solution, at a time t > to, after a singularity has 
appeared. The conservative solution has two strong singularities at 7“(t) < 7''’(t), while the dissipative 
solution has a strong singularity at 7“(t) and a weak singularity at 7(t). The two solutions coincide 
for X < 7“(t) and for x > ^{t). 
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2. For X > Xq we denote by y = a'^{X) the curve where re = vr, in the dissipative solution. 
A Taylor approximation for is derived from the identities 

w{X,Yo) = wo + wxx{Xo,Yo) • ~ + 0(1) • (A - Ao)^ 

wy{X,Y) = w(Xo,yo) + 0(l)-(|^-^o| + |h^-i"o|), 
valid in the region where w < tt. Together, they imply 

u#(A) = Yo + K{X-Xof + 0{l)-{X-Xof, (4.6) 

where k > 0 is the same constant found in (3.37) for the conservative solution. 

For Y' > Yq, (3.38) and (4.6) together imply 

/w/ _ Y \ 

X^{Y')-X-{Y') = 2 I -+0(1) • ly'-Tol- (4.7) 



Figure 8 ; Estimating the values of a dissipative solution near a singularity. Notice that the functions 
X, t, u are constant on every horizontal segment contained in the shaded region where w = ir. 

3 . Consider a point {X,Y) with X > Xq and Y < a^{X). By the second equation in (2.9) it 
follows 

uiX, Y) = uiX, Yo) + £ (?) (A, y') dY'. (4.8) 

As in Fig. 8, for Y' G [yo,y], call X~{Y') and A**(y') respectively the points where a~{X) = 
Y' and cj^(A) = Y' . Since zx = Qx = ^ when w = n, hy the second equations in (4.1) and in 
(4.2) we have 

ziX,Y') = ziX-iY'),Y')+ r zx{X',Y')dX' 

Jx-{Y') 

= ziX-iY'),Y')+ f (cos w-cos z)p) {X',Y')dX\ 

JxHY') \8c^(u) ) 

q{X,Y') = (?(A-(yO,y')+ r (4^ [smw-smz]pq) {X',Y')dX'. 

JxHY') \8c^(u) ) 

4. For notational convenience, in the following we denote by {x,t,u,w, z,p,q){X,Y) the 
components describing a dissipative solution, and by (x, t, u, w, z,p, q){X, Y) the corresponding 


(4.9) 

(4.10) 
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components of the conservative solution. We observe that all these functions are Lipschitz 
continuous. As shown in Fig. 6, these two solutions can be different only at points {X,Y) in 
the region bounded by the curves a~ and a, namely 

{X<Xo, Y>a-{X)} U {X>Xo, Y>Yo}. 

Consider a point {X, Y') with X > Aq, Y' < a^{X). By (4.1), observing that z = z for Y < Yq 
and using (4.7) we find 

z(X,Y')-z(X,Y') 


/ zx{X', Y') dX' + {zx- zx) (A', Y') dX' 

Jx-{Y') Jxi(Y') 


X 
Xtt(Y') 


+0(1) • (A«(y') - A-(y'))^ + 0(1) • {X - A“(y')) (y' - yo) 


_ c'(iio)(l + COSZo) 
4c^{uo) 

By (4.2), a similar computation yields 
q{X,Y')-qiX,Y') 


Po 


• (y' - Yo)Y^ + 0(1) • |y' - Yo 


(4.11) 


rXO(Y') j-X 

/ qx{X',Y')dX'+ / {qx-qx){X'X)dX' 

lx-(Y') Jxt{Y') 




(4.12) 


+0(1) • (A»(y') - A-(y'))^ + 0(1) • (A - A»(y')) (y' - yo) 


c'{uo) sin zq 
4c^(uo) 


poqo ■ (y' 


yo)l/2 + 0(1) . |y/ _ 


Next, using the second equation in (2.9) and recalling that u{X,Yo) = u{X,Yo), for any 
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X G [Xq, Xq + 1] and Y G [Yq, a^{X)] we obtain 


u(x,r)-u(x,r) = 

- £ (i5 - " SS ■« -> - ^ <“ - “0 <"■ 


+0(1) ■ j {\z- z|2 + \q - q\^ + \u- uf^ (X, Y') dY' 

+O{l)-(^\X-Xo\ + \Y-Yo\y £ y-z\ + \q-q\ + \u-u\)iX,Y')dY' 

= rjo-{Y-Yof/^ + 0{l)-\Y-Yo\^ 
where the constant r/o is computed by 


(4.13) 


2 r c'(tto)(l + cos zq) 1 coszo 2 \ c'{uq) sinzo 1 sinzo 

5. Using the second equations in (2.12) and in (2.13), we obtain similar estimates for the 
variables x,t. Namely, 


x{X,Y)-x{X,Y) = - 


1 + cos z 1 + cos z 
-^- Q -^- 


q)iX,Y>)dY> 


sinzo .. . 1+coszo \\ r v vf\ jv> 

qo-{z-z) - - - ■ {q - q) j {X,Y ) dY 


+ 0 ( 1 ) ■ £ y- + Iq - (^> ^0 dY' 

+0(1) • (|X - Xol + |y - Uol) • £ {\z - z| + \q - g|) (X, +') dY' 


(4.15) 


= O(l)-|F-Fo 
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Indeed, the coefficient of the leading order term 0(1) • {Y — vanishes. Similarly, 

t-Y 


t{X,Y)-t{X,Y) = - 
_ f sin zo 


'Yo 


1 + cosi^ l + cosz 


4c(u) 


■q 


4c(ii) 


q]iX,Y')dY' 


qo-{z-z)- ^ I -{(1-0)+ • (« - ^^) ) {X, Y') dY' 


Iyq \4c(rto) 4c(uo) 4c2(rio) 

X 

'Yo 

rY 
'Yo 


+0(1) • (|l - z\^ + \q- q\^ + 111 - u|2) (X, y') dY' 

+0(1) • (|X - Xol + |y - yol) • £ [\z -z\ + \q-q\ + \u- u|) (X, X') dY' 


= o(i)-|y-yc 


01 


(4.16) 


6. The estimate (4.13) provides a bound on the difference ii — u between a conservative 
and a dissipative solution, at a given point (X, X). However, our main goal is to estimate 
the difference u — u as functions of the original variables x,t. For this purpose, consider a 
dissipative solution u and a point 

p = (x,t) = (x(x,y),t(x,y)), (4.17) 

with 

X > Xo, Xo < y < a\x). (4.18) 

Moreover, let u be the conservative solution with the same initial data, and let (X,X) be the 
point which is mapped to P in the conservative solution, so that 

P = (x,t) = (x(X, X),t(X,X)). (4.19) 

Using (4.13), (4.15), (4.16), and recalling that the conservative solution u = u{x,t) is Holder 
continuous of exponent 1/2 w.r.t. both variables x,t, we obtain 

\u{x,t)-u{x,t)\ < |u(X,X)-n(X,X)| + |u(X,X)-u(X,X)| 


= 0(1) • (|x(X,X) - x(X,X)|'/' + |t(X,X) - t(X,X)|'/') + 0(1) . |X - Xo|3/2 
= 0(1) • (|x(X,X) - x(X,X)|'/' + |t(X,X) - t(X,X)|'/') + 0(1) • |X - Xo|3/2 


= 0(l)-|X-Xo 


In a neighborhood of (xo,to) we have 




1 + cos z 1 + cos Zq 


(4.20) 

(4.21) 


For (x,t) as in (4.17)-(4.18), one has 

t-to = [f(X,X)-t(X,Xo)] + [t(X,Xo)-t(Xo,Xo)] > ^^^^'io-lY-Yo 
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Together with (4.20), this proves (4.5). 


7. It remains to prove that the solution u = u{t, x) is not twice differentiable along the forward 
characteristic 7 . 

Consider a point (xi,ti) = [x{X,YQ),t{X,YQ)) on 7 , with X > Xq. Let x = ^{t) be the 
backward characteristic through (xi,ti), so that 


7 (ti) = xi, 7(i) = c(tx( 7 (t),t)). 


Assume that u were twice differentiable at the point (xi,ti). Then the map t 1 -^ u{'y{t),t) 
would also be twice differentiable at t = ti. Indeed 


dt 


= -c{u)ux + ut. 


= -c'{u){-c{u)ux + ut)ux + c^iu)uxx-‘2c{u)uxt + utt- 
To reach a contradiction, consider the map t Y (r) implicitly defined by 

t{X,Y{T)) = T. 

By (4.21) this map is well defined. In particular, Y(ti) = Yq. We thus have 

4c(ii) sinz 


(4.22) 


= ^u{X,Y{t)) = UY 


dt 


(l + cosz )(7 1 + cosz 


(4.23) 


We now show that this first derivative cannot be a Lipschitz continuous function of time, for 
f ~ ti. Indeed, by (4.23) and the mean value theorem we have 


Ux,Yit))-Ux,Y(h)) - 

dt dt 1 + cos z(A, y(t)) 1 + cosz(a, To) 

1 


(4.24) 


1 + cos z 


[z(A,T(f))-z(X,To)], 


for some intermediate value Y^ G [Lo,T(t)]. Call z = z(A, T) the corresponding conservative 
solution. Observe that z is smooth and coincides with z on the horizontal line {Y = Tq}. 
Using (4.11) we obtain 

|z(X,T(t))-z(X,To)| > [z{X,Y{t))-z{X,Y{t))]-[z{X,Y{t))-z{X,Yo)] 

> ^ (y(t) -- 0(1) ^ |y(i) - ni. 

As a consequence, for t^ti, the function t i-)- z{X,Y{t)) is not Lipschitz continuous, and the 
same applies to the left hand side of (4.24). We thus conclude that the map t i-)- u{'y{t),t) 
cannot be twice differentiable at f = U, in contradiction with (4.22). This completes the proof 
of Theorem 4. □ 
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